Role of Diffusion in Two-dimensional Bimolecular Recombination 
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Experiments on carrier recombination in two-dimensional organic structures are often interpreted 
in the frame of the Langevin model with taking into account only the drift of the charge carriers 
in their mutual electric field. While this approach is well justified for three-dimensional systems, 
it is in general not valid for two-dimensional structures, where the contribution of diffusion can 
play a dominant role. We study the two-dimensional Langevin recombination theoretically and 
find the critical concentration below which diffusion cannot be neglected. For typical experimental 
conditions, neglecting the diffusion leads to an underestimation of the recombination rate by several 
times. 



The problem of charge carrier recombination is one of 
the central problems in the solid state physics since re- 
combination processes determine all optoelectronic prop- 
erties of a device or a material under study. For in- 
stance, recombination can be an important factor lim- 
iting the conversion efficiency in organic solar cells. In 
low-mobility organic semiconductors charge carrier re- 
combination is usually dominated by the Langevin bi- 
molecular mode. 1 Langevin recombination is character- 
ized by the rate limiting step being finding the opposite 
charge carrier, independent of the subsequent recombi- 
nation mechanism. It is well-known that in the three- 
dimensional (3D) case the Langevin recombination rate 
is determined by drift motion of charge carriers due to 
their mutual Coulomb attraction and that diffusion can 
be neglected. 2 The recombination rate per unit volume 
equals (fi e + ^ft)eA r2 /(ee ), where N e is the electron den- 
sity which is assumed to be equal to the hole density 
A e = Nh = A, n e and fih are the electron and hole 
mobilities, respectively, esq is the dielectric permittivity, 
and e is the electronic charged— This result is obtained 
by integrating the drift flux of electrons through a spher- 
ical surface of radius r surrounding a hole. Since the 
electric field scales as r~ 2 and the surface area of the 
sphere scales as r 2 , the value of r chosen is unimpor- 
tant, leading to a simple solution with constant electron 
density, thus justifying the neglect of diffusion^ The re- 
combination rate per unit volume is often written in the 
form /3A 2 via the bimolecular recombination coefficient 
ft = e(/i e + /x^)/££o, independent of carrier concentration 
in 3D. 

Numerous recent experimental studies on bulk hetero- 
junction solar cells have demonstrated that the recom- 
bination in organic materials with lamellar structures 
is much slower than what is predicted by Langevin's 
theory^— Hence, in these films, where 2D conductivity 
mode is expected, the bimolecular recombination coef- 
ficient depends on the charge carrier density N£&H-~— 
The straightforward reformulation of the Langevin 3D 
formalism for the 2D case indeed gives a concentration- 



dependent recombination coefficient /?2D oc A 1 / 2 if only 
the drift motion is taken into account. This result was 
considered as an explanation of the experimental data^^ 

It is easy to understand the latter result rewriting 
the recombination rate per unit volume in the form 
TZ-rec = N/t, where t is the so-called recombination time. 
Remarkably in any model that neglects diffusion, the re- 
combination time in 2D is proportional to A -3 / 2 , and 
consequently the recombination rate per unit volume is 
proportional to A 5 / 2 . It follows from the fact that there 
is only one combination of the parameters e, [i and A 
(not related to diffusion), which has the dimensionality 
of time: (e e/efi)N~ 3 / 2 . 

It is however not possible to neglect the contribution 
of diffusion to the bimolecular Langevin recombination 
in 2D systems^ In two dimensions, the electric field still 
scales as r -2 , but the integration takes place over a circle 
of circumference 27rr. To achieve a recombination current 
which is independent on r requires a density of electrons 
which varies with r. The presence of a carrier density 
gradient means that diffusion must be explicitly included 
in the problem^ 

The combined effects of the drift and diffusion on the 
bimolecular recombination in 2D has been studied by 
Greenham and Bobberti^ In the following we first briefly 
repeat their arguments. Afterwards we solve the ob- 
tained equations analytically and compare the results 
with Ref. ! Finally we reformulate the model of Ref. 
in order to make it compatible with the model considered 
in Refs. d|3 in order to estimate the effect of diffusion. It 
will be shown that the effect depends essentially on the 
total concentration of charge carriers. 

The following model was considered by Greenham and 
Bobbertj- 2 Electrons come uniformly to a circular area 
around a hole. The radius R of the circle is chosen so 
that its inverse area is equal to the hole concentration 
N h : 



(1) 



The hole at the origin serves as a drain for incoming 



2 



electrons. Electrons can move within this area, due to 
drift in the electric field of the hole, 



E(r) 



47T£ £ r 3 



(2) 



and also due to diffusion. 

This model can be formalized via the following set 
of equations for the electron flux density J(r) and the 
coordinate-dependent electron concentration n(r): 

V • J(r) = /, (3) 
J(r) = -Am(r)E(r)--DVn(r), (4) 

that is fulfilled in the range < r < R. Here D is 
the diffusion coefficient (D = fikT/e according to Ein- 
stein's relation), / is the incoming flux of electrons. The 
boundary conditions for this problem are the following: 
J(R) = (there is no electron flux out of the circular 
area), and n(0) = (the hole is an ideal absorber of elec- 
trons). Solving this set of equations, one can obtain the 
recombination time t as 



t = (n(r))/f, 



(5) 



where (n(r)) denotes the concentration averaged over the 
circular area. Then one can obtain the recombination 
rate 7^-7*g£ as 



K rec = N e /t, 



(0) 



where N e is the total electron concentration in the sys- 
tem. 

The result of Greenham and Bobbert£ can be conve- 
niently represented in terms of the dimensionless param- 
eters 



a 



D 



where a = e 2 / ATTe^ekT is the Onsager radius: 
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where G is the Meijer G function. From equations ((T|), 
©, ([7]), and ([5]) one can obtain the concentration depen- 
dence of the recombination time. 

Greenham and Bobbert^ analyzed the obtained result 
numerically and suggested for the concentration depen- 
dence of the recombination time a power-law interpola- 
tion: t ~ iV^ 1 ' 43 . (As mentioned above, we will assume 
that electron and hole concentrations are equal: N e = 
Nh = N.) Using Eq. © one obtains then the following 
dependence for the recombination rate: lZ rec ~ N 2A3 . 

One should however mention that the above interpre- 
tation of the numerical result is valid only in a restricted 
range of concentrations. In order to obtain a more gen- 
eral picture, let us consider the limiting cases of large 



and small concentrations. For this purpose let us use the 
following series expansions of the Meijer G function: 
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+— +0{x 4 ) at x -+ 0, (9) 
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o(l) at x — >• 00. 

(10) 

Here 7 sa 0.577 is Euler's constant. 

Let us first consider the limit of large concentrations, 
i. e. small R. Taking into account the expansion ([§]), it 
is easy to show that 



t=^R 3 + 0(R 4 ), 
whence the dependence t(N) follows the relation: 



t = 



15n 3 / 2 aD 



N -3/2 = 



£0£ ^-3/2^ 
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(11) 



(12) 



Consequently, lZ rec ~ N 2 - 5 , in accordance to the result 
of Juska et a/.,— where only drift was taken into account 
and diffusion was neglected. Remarkably, in the case of 
high carrier concentrations N the contribution of diffu- 
sion to the recombination rate is not important. Notably 
our coefficient of proportionality (Ib^/TieiJL/SeQe) appears 
2.5 times larger than the one given in Ref. |3j. This differ- 
ence comes from different assumptions in the models: in 
Ref. [H each electron was considered to start at the same 
distance R from the hole, while in the model considered 
above (Ref. Q), the starting points of electrons were cho- 
sen randomly within the circle of radius R. 

We reformulate the formalism of Greenham and 
Bobbert2> adjusting it to the model assumptions of Juska 
et a/M To do so one has to replace the uniform influx of 
electrons inside the circular area by the influx from only 
the periphery of the circle. Then Eq. Q takes the form: 



V-J(r) =0, 



(13) 



with replacing the boundary condition J(R) = by the 
relation J Z (R) = —fR/2. This condition provides the 
same average influx of electrons as Eq. ([3]). The solu- 
tion of this modified problem gives instead of Eq. ((SJ) the 
following expression: 
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In the limit of small R (i. e. high concentrations) one can 
keep only the first term of the expansion © and obtain: 



t=~R 3 + 0(R% 



and consequently 
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FIG. 1: The ratio between the recombination times thrift ob- 
tained by taking only drift into account [Eq, (ll6|l ] and idrift+dia 
obtained when considering both drift and diffusion [Eq. (|14|l ]. 
as a function of the charge carrier concentration N. The con- 
centration is measured in units of the number of charge car- 
riers inside the Onsager radius a. The arrow corresponds to 
the experimental concentration in Ref. 0. 



which is exactly the result of Juska et al£& Figure [T] illus- 
trates that at high carrier concentrations one can indeed 
neglect the contribution of diffusion to the recombination 
time. However, the smaller is the concentration the more 
important is the contribution of the diffusion process. For 
instance, at carrier concentrations in the experiments re- 
ported in Ref. neglecting the contribution of diffusion 
gives a recombination time three times longer than when 
diffusion is included. At lower concentrations this over- 
estimation can reach orders of magnitude. Fig. [T] shows 
that within the accuracy of 50 % the neglect of diffusion 
is justified for concentrations TV > N c = In other 



words, if the average distance between an electron and a 
hole is larger than the Onsager radius a, which is the case 
at low charge carrier concentrations, one cannot neglect 
the contribution of diffusion to the Langevin recombina- 
tion in 2D systems. 

Let us therefore pay more attention to the limit of low 
concentrations. Using equations (JTJ, ||7J), ©, (ITU1) . and 
(fTH) one comes to the conclusion that in this case the 
concentration dependence of the recombination time has 
the form: 



' TV- 1 log " 



4nD 



a 2 A' 



(17) 



where the numeric constant a is l/(7re 27+3 / 2 ) ps 0.022 for 
uniform generation of carriers inside the circular 
considered in Ref. 0, and l/(7re 27+1 ) « 0.037 for genera- 
tion of carriers on the periphery of the circle, as consid- 
ered in Refs. In the limit N — > 0, the recombination 
rate reads 



_ N 



AttD N 2 flog ■ 



a 



t . V-°a>NJ. ■ (18) 

Using for the recombination rate per unit volume the 

traditional notation (3N 2 we conclude that in the limit 
of low concentrations N, the bimolecular recombination 
coefficient (3 depends logarithmically on N. This depen- 
dency is essentially weaker than the one predicted by 
considering solely the drift processes.^ 
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